Berry Phase of Nonlinear Correction 
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We investigate the geometric phase or Berry phase of adiabatic quantum evolution in the Bose- 
Einstein condensate (BEC) systems governed by nonlinear Gross-Pitaevskii(GP) equations. We 
study how this phase is modified by the nonlinearity and find that the Bogoliubov fluctuations 
around the eigenstates are accumulated during the nonlinear adiabatic evolution and contribute a 
finite phase of geometric nature. A two-mode BEC model is used to illustrate our theory. Our 
theory is applicable to other nonlinear systems such as paraxial wave equation for nonlinear optics 
and Ginzburg-Landau equations for complex order parameters in condensed-matter physics. 
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Adiabatic theory, as a fundamental issue of quantum 
mechanics, involves two aspects: i) when the Hamilto- 
nian changes slowly compared to the level spacings, an 
initial nondegenerate eigenstate remains to be an instan- 
taneous eigenstate[l|; ii) the phase acquired by the eigen- 
state is the sum of the time integral of the eigenenergy 
(dynamical phase) and a quantity independent of the 
time duration and related to the geometric property of 
the closed path in parameter space (geometric phase or 
Berry phase) 0, • The adiabatic theory has been play- 
ing a crucial role in preparation and control of quantum 
states IH • Recently, the Berry phase and related geomet- 
ric phasesf^, Q has received renewed interest due to its 
important use in implementation of quantum computing 
gates and applications in condense matter physics Q. 

For a nonlinear quantum system, such as that de- 
scribed by the nonlinear Schrodinger equations, how the 
adiabatic theory gets modified. Nonlinear quantum sys- 
tems have become increasingly important in physics. 
They often arise in the mean field treatment of many- 
body quantum systems, such as Bose-Einstein conden- 
sates (BECs) of dilute atomic gases 0. Recently, extend- 
ing the first aspect of adiabatic theory to the nonlin- 
ear systems, i.e., investigating the adiabatic condition 
and adiabaticit y f or the nonlinear quantum evolution 
has been done |ldl . Hlj ]. interestingly it was found that 
the adiabaticity of an eigenstate only requires that the 
control parameters vary slowly with respect to the Bo- 
goliubov excitation frequencies and has nothing to do 
with the level spacings. Nevertheless, Berry phase is- 
sue in such nonlinear system is far from well understand, 
while some superficial observations suggest that Berry's 
formula without any correction is applicable to nonlin- 
ear system suppose the system is invariant under gauge 
transformation of the first kind [13] • 

In this letter, we make a thorough analysis on geomet- 
ric phase associated with the adiabatic evolution of an 
eigenstate, and strikingly we find that the Berry phase is 
dramatically modified by the nonlinearity. This finding is 
completely contrary to previous superficial observations. 
The underlying mechanism has been revealed: for a non- 



linear system because the Hamiltonian is a functional of 
the instantaneous wavefunctions, the Bogoliubov fluctu- 
ations around the eigenstate caused by the slow change 
of the system are allowed to feedback to the Hamiltonian. 
They are accumulated during an adiabatic evolution and 
eventually contribute a finite phase of geometric nature. 
A two-mode BEC model is used to illustrate our theory. 

Let us consider the BEC trapped in a potential 
V(R; r), the evolution of its wavefunction is governed by 
following nonlinear Gross-Pitaevskii(GP) equation (h = 
m = 



(i) 



with Ho = —4 V +V(R;r), where R is the parame- 
ter vector, g is the nonlinear parameter representing the 
interaction between the coherent atoms. The total en- 
ergy of the system Er = J drE(ip* ,ip), where the en- 
ergy density E(tf>,ij)*) = tp*Hotp + ^g\ip\ 4 . The above sys- 
tem is invariant under gauge transforms of the first kind, 
ip(r, t) —>■ exp(i7y)-0(r, t) with constant rj. The gauge sym- 
metry implies that the total atom number is conserved, 
i.e., ] dr\yb\ 2 = 1. 

Let A be the overall phase of the wave function; we 
may take it to be the phase of the wavefunction at a 
fixed position tq, for example, A = — arg(^(ro, t)). We 
split off this overall phase by writing ip = e~ lX cf>, then cf> 
belongs to the so called projective Hilbert space. From 
([]]) we obtain 

' /A ' "■ ■ f M(f,^S#'#). (2) 



dt 



The eigenequation of the system is 
H tp + g\ilj\ 2 tp = if>P- 



(3) 

where ip is the eigenfunction and \x is the eigenvalue (or 
chemical potential). 

Now we consider the parameter vector R varies slowly 
in time, and introduce the dimensionless adiabatic pa- 
rameter of e ~ |^| as the measure how slow the param- 
eters change. The adiabatic parameter tends to zero, i.e., 
£ — * 0, indicating the adiabatic limit. 
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Consequently, the expression of the total phase can 
be expanded in a perturbation series in the adiabatic 
parameter, i.e., 



^ = a (e°) + ai(e 1 ) +o(e 2 ). 



(4) 



When the parameters move in a circuit, the eigenstates 
evolves for an infinite long time duration in the adiabatic 
limit. The time integral of the zero-order term gives so- 
called dynamic phase because it is closely related to the 
temporal process of the evolution. The time integral of 
the first-order term gives an additional contribution to 
the overall phase, which will be shown later is of a geo- 
metric nature, that is, only depends on the geometry of 
the close path in the parameter space. The contribution 
of the higher-order term vanishes in the adiabatic limit. 

In the quantum evolution with slowly-changing param- 
eters, we assume <fi — <^(R) + S(f>(R), where 4>(R) is the 
wavefunction of the instantaneous eigenstate correspond- 
ing to the local minimum energy. 8<j) denotes the secu- 
lar part of the Bogoliubov fluctuations induced by the 
system's slow change while the rapid oscillations in the 
fluctuations are ignored because they vanish after a long- 
term average. 8<p depends on the adiabatic parameter 
and is of order s, then from Eq.([2]) and with the help of 
relation Eq.(3), we have the explicit expressions as fol- 
lows, 



a (e°) =/x(R), 



(5) 



ai(e 



ill 



5<(0 2 r, ^)\{84>M*)), (6) 



where we denote (a\b) = J dr{a*b),({a,b)\{c,d)) = 
J dr{a*c + 6*<i),and the second order term like ^^has 
been ignored. 

From the above expressions we see that the dynamical 
phase has been modified to be the time integral of the 
chemical potential rather than the energy. This is be- 
cause the instantaneous eigenstates are feedback to the 
Hamiltonian. More interestingly, the first-order term, 
i.e., the Berry phase term has been modified due to the 
feedback of the Bogoliubov fluctuations to the Hamilto- 
nian. To evaluate it qualitatively and express the modi- 
fied geometric phase explicitly, let us introduce a set of 
orthogonal basis \k) and the variable ipj is the j-th com- 
ponent, i.e., ipj — (j\ip). Without losing generality, the 
projective Hilbert space is set to be of a specific gauge 
that the phase of the N-th component is zero. In the 
projective Hilbert space, the new variables (rij,9j) are 
introduced through <frj = ^/nje iej . Substituting the ex- 
pression of ipj = y/nje 161 ' e~ l S ® dt into the GP equation, 
and separating real and imaginary parts, we have follow- 
ing differential equations for the density rij and phase 9j , 
respectively, 



drij „ d9j 

—r- = fi,-r-= j = 1,2,. 
dt JJ ' dt Jl J 



.N — 1 



(7) 



where 



N 

fe=l 



Q fe (R)e i(e 



+2g 



N 

E 

k.l.m— 1 



(8) 



hj = -y J— Re C jk (R)e l 

k=l 



N 



> E 

k.l.m— 1 



■P, (9) 



P 



N 

E 

k=l 



Re 



c Nk (Ry { 



e k -e N ) 



N 



^ E 



k.l.m— 1 



n N 



](10) 



with C jk (R) = {j\H {B.)\k),D jikt i, m = (j\(k\l)\m). 

The last equation is from the fact that in the pro- 
jective Hilbert space the phase of wavefunction of the 
N-th component is set to be zero. The norm conver- 
sation condition = 1 — Y^k=i n k could be used to 
remove the variable n^r in the above equations. In the 
representation of new variables, (jij,9j) satisfy equations 

of the equilibrium state, i.e., f^f-, ^jjM |^ g^ = 0. 

(rij,9j) are function of the parameter R corresponding 
to the eigenstates of GP equation. Let us make per- 
turbation expansion around the eigenstate with rij — 
rij(R) + Srij(R),9j = 0,(R) + S9j(R). Here, ^-(R) = 

v /n i (R)e i ^W, (6nj(K),S9j(K)) are the fluctuations de- 
pending on the adiabatic parameter and of order e. Then 
inserting the above expansion into the equations (7) and 
ignoring the higher order terms such as , , with 
denoting v — (m, 9\\ n^-x, 0/v-i) wc obtain that 

dv dR 
dR~dt 

where the matrix takes the form 



CSv, 



(11) 



£ = { £ 3k} {N _ i _ N _ i) ,£ 3 k = 



dn k d6 k 
dhj dhj 
dn k d6 k 



then, inversely we have, 



5v = C~ l 



dv dR 



(12) 



(13) 



dR dt 

The differential relation between the new variables and 
old ones take the form 

Sri* 



o - 



n, 



(14) 



3 



in which 



/ In- 1/2 e^ 

n - I 2 j , 

3 ~ 1— -1/2 -id- ■ I— -ih 



(15) 



Substituting (fTB"]) and (TTJI into © , we finally obtain the 
explicit expression of adiabatic geometric phase that con- 
tains two terms, 

lg = lB+ 1NL (16) 

where the first term is the usual Berry phase formula, 



IB 



N—X — 

3=1 



and additional term is from the nonlinearity, taking the 
form, 



Inl =9 j> ^A|IIo£ 



dR. 



(18) 



Here, A 

N-l 

E % - 

3=1 



(fa 



E «j 

3=1 



1)R 



IK 



(njv- 



_l/2 0, ,_ 



AT-1 

E % 

3=1 



l)n 



1/2 
JV— ] 



(n N - 



dv / dni d6 

da ~ y i 
n 



N-l 

- E «3 
3=1 

dnj\r-i d(?jv-l\T 



in— 1/2 _ 



, and diagonal matrix 



dR ' dR ' dR ' dR 

diag(IIi , II2 , . ..n^v—i) -Notice that to simplify the 
expression of A, we use the approximation that the 
overlap intergral Dj^j,™, — when the subscripts are 
not all identical. 

Both 7s and jnl have the geometric property of pa- 
rameter space. The novel second term indicates that, 
the Bogoliubov fluctuations induced by the slow change 
of the system that is negligible in linear case, however, 
could be accumulated in the nonlinear adiabatic evolu- 
tion and contribute the finite phase of a geometric na- 
ture. This is because the Hamiltonian in the nonlinear 
system contains the instantaneous wavefunction and as a 
result the Bogoliubov fluctuations around the eigenstate 
are allowed to feedback to the Hamiltonian. 

As an illustration of our theoretical formulism, we con- 
sider the two-mode BEC model as example, which is de- 
scribed by the following GP equation [14j, 



with 



dt V * 



ff(*i,*s 



*1 
*2 



Z\*!\' 

Ppitp 
2 e 



P p -iip 



(19) 



(20) 



the R = (p, (p, Z) are parameters. For simplicity, we 
fix p and Z and change the parameter <p from to 2ir 
adiabatically. 



The energy of the system £ , ($i,$ 2 ) = f (|*i| 4 + 
|*2| 4 ) + § (e-^W^a + e^*i*J)-The above system is 
invariant under gauge transforms of the first kind. 

£1 

$ 2 



The eigenequations take the forms of H 



For the nonlinear system, the number 



' i>2 

of the eigenstates may be larger than the dimen- 
sion of the Hilbert space and the eigenstate could be 
unstable [nj. We have obtained four eigenstates for the 
case Z > p by solving the above eigenequation. Three 
of them are stable and one is unstable. We choose 

following stable eigenstate for example to illustrate 

1/2 

i(l-Jl 



our theory, i.e.,*&i 



£1 



1 



1 



1/2 



e %v , with eigenvalue of p = Z. 



We choose the overall phase A 
equation (fH)]) , we have 



dX 
~~dl 



d 



arg (^1) , from GP 



Z 



*((*!, $2)|7d($l,$2))+£($l,$ 2 ) + ^(|$l| 4 +|$2| 4 ). 



'dt 



(21) 

In the adiabatic evolution, we assume $i = &i{(f) + 
5$>i((p), where <5$i is the fluctuation depending on the 
adiabatic parameter of order e = dip/dt. Then we have 
the explicit expression of the zeroth order, 



a (e°) = Z, 
and using conservation of the particle |$i| 



(22) 



0,we have the first-order term, 



l,and the fact arg$ x 

+Z((2|¥ 2 | 2 - 1)$;,(2|$ 2 | 2 - l)$ 2 |(<y$ 2 ,$*;))(23) 



at 



Berry term of the geometric phase is readily deduced, 

lB=n ^l+^l-^y (24) 

Now we are going to derive the additional term which 
gives the jnl of the geometric phase. Let us in- 
troduce the new variables (n, 0) through ($i, < I ) 2 ) = 
(VT"^, y^e 19 ). Substitute (*i,* 2 ) = e^f* 0dt ($ u $ 2 ) 
into Eq. (|19p . and separate the real and imaginary parts, 
we have four differential equations, in which two of them 
are identical due to the norm conversation, 



dt 
d 
~dt 



—py n — n 2 sin(# — ip) 
p\J\ — n 



2V^ 



ip) — Zn + [3, 



P = 



Z{\-n) + Uj^ux{9-<p). 



(25) 
(26) 
(27) 
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The eigenstate make up the fixed point of equations (|25|) 
and flU, i.e., n = 3(1 + \/l-£),6 = <p. Let us 



make perturbation expansion around the eigenstate with 
71 = ri(ip) + Sn,9 — 9(ip) + 89. Then inserting the above 
expansion into the equations (j25jl and (j26jl . and ignoring 



(28) 



the higher order terms of we obtain that 




with 




(29) 



1NL 



f( A|n£ " 



dn 89 

dip' dtp 



dip. 



(30) 



where A = Z((2n— \)\Jn , (2n — l)yn ). Finally, we 
get 



7rp 



(31) 



The above new form of geometric phase is verified nu- 
merically by directly integrating the Schrodinger Eq.(19) 
using Runge-Kutta algorithm. 

In summary, we show that the Berry phase associated 
with the adiabatic evolution of BECs has been modi- 
fied by the nonlinearity. This nonlinear correction is of 
significance because it could affect the interferences of 
the matter waves. It is expected to be observed in fu- 
ture's experiments. Our theories, on the other aspect, 
cover the general nonlinear systems in nonlinear optics 
and condense matter physics. Since some of the theories 
covered by our results are meanfield limits of quantum 
many-body theories, the possibility of generalizing these 
considerations to quantized field theories from the corre- 
spondence principle is of great interest for future study. 
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